Wigner time-delay distribution in chaotic cavities and freezing transition 
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Using the joint distribution for proper time-delays of a chaotic cavity derived by Brouwer, Frahm 
& Beenakker [Phys. Rev. Lett. 78, 4737 (1997)], we obtain, in the limit of large number of channels 
N, the large deviation function for the distribution of the Wigner time-delay (the sum of proper 
times) by a Coulomb gas method. We show that the existence of a power law tail originates from 
narrow resonance contributions, related to a (second order) freezing transition in the Coulomb gas. 

PACS numbers: 05.60.Gg ; 03.65.Nk ; 05.45.Mt 



The study of scattering theory in chaotic or disor- 
dered quantum systems within the random matrix theory 
(RMT) has been a subject of intense research for many 
years. Though originated in nuclear physics (see the re- 
view [I]), it has major implications in condensed matter 
theory where it can be used to model electronic transport 
in mesoscopic (coherent) conductors The dynam- 

ics of an electron of energy E is described through the 
N x N on-shcll scattering matrix S(E), where N is the 
number of scattering channels. A useful concept that 
characterises the temporal aspects of the scattering pro- 
cess is time-delay [U [5] undergone by an incident wave 
packet. This is captured by the Wigner-Smith time-delay 
matrix 6 1, Q{E) = -iS{E)^ (with % = 1), whose 

eigenvalues are the proper time-delays t\, ■ ■ ■ , Ty. 

If the system is characterised by some complex dy- 
namics, due to the presence of disorder or chaos, the sta- 
tistical properties of time-delays exhibit interesting uni- 
versal characteristics : the universality of the time-delay 
distribution for ID-disordered quantum mechanics was 
demonstrated in [7] (see also [5HT2], and [J3] for 2D & 3D 
cases). The situation where the dynamics is chaotic has 
been extensively studied within RMT: the marginal law 
of partial time-delays [T3], §n{t) — -h J2 a (^( r — T a )), 
was obtained for GUE symmetry (Dyson index /3 = 2) in 
[Tfjl HB] . The time-delay distribution was obtained in the 
N = 1 case with G. {1. 2. 4} in [2"0]. Using the "alter- 
native RMT" introduced in [21j . Brouwer and coworkers 
succeeded in finding the joint distribution of the inverse 
proper time-delays 7fe = 1/Yfc [HI [33] : 
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(the times are measured in units of the Heisenberg time 
th = 2ttH/A, where A is the mean level spacing). This 
measure, known as the Laguerre ensemble of random ma- 
trices, also corresponds to the distribution of the (posi- 
tive) eigenvalues of Wishart matrices X'X, where the 
matrix X has size N x (27V — 1 + 2//3) with i.i.d. Gaus- 
sian matrix elements. 

In this article we are interested in the Wigner time de- 
lay, defined as the sum of proper (or partial) [T3] time 



delays r w = f ±Tr{Q} = ^^2 a=1 T a . This quantity is 
of great interest due to its close relation to the density 
of states (DoS) of the open system, through the Krein- 
Friedel relation [T7] : v(E) = ± Tr{Q(£:)} = ^Nt w . 
The Wigner time delay (or related quantities such as in- 
jectance or emittance) is a central concept for studying 
charging effects, e.g. for mesoscopic capacitances [HI [20]. 
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FIG. 1: (color online). Sketch of the distribution of s = Ntw- 
The dashed line at s = sm — 1 + (4/iV) 1,/3 separates the two 
phases of the Coulomb gas with densities represented in the 
small figures on the left and the right respectively. 



We denote by Pjv(t) = (6(t - ^ J2 a T i)> the Wigner 
time-delay distribution. Despite the fact that the joint 
distribution of proper times is known already for 15 years, 
little is known about the distribution of rw for general 
N: it has been computed explicitly only for N = 1, 
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3-57 [20] and N = 2, P 2 (t) = 
where U(a, b; z) is the confluent hypergeometric 



The distribution was conjectured to have a 
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power law tail for large r, Pjv(t) ~ t " 2 " in _ 
(for j3 = 2) by using the resonance picture allowing to 
identify the tails of Pn{t) and Pn(t) (for a heuristic 
argument using the relation to resonance width, cf. the 
review [25]). More recently, the first three cumulants 
of T\y were derived by a generating function method 
[2"5] . However a full understanding of its distribution for 
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general N is still missing so far. 

In this Letter, by analysing an underlying Coulomb gas 
we provide a complete description of Pn(t) for large N 
and show that it has a rather rich behaviour including 
an interesting nonanalytic point which is a consequence 
of a freezing transition in the Coulomb gas. Limiting 
behaviours of Pn(t) may be summarised as follows (t\v 
is measured in unit of th) : 
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A sketch of the distribution is given in Fig. [TJ The Gaus- 
sian form around r ~ 1/N in ([3| allows one to extract 
the mean time-delay and its variance. Reinstating th, we 
obtain (rw) = ^y- Consequently, the mean DoS reads 
{v(E)) — N{tw)/2tt — 1/A as expected. Similarly, the 
variance can be read off ([3| 
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Eq. @ agrees with the leading order of the result ob- 



tained in Ref. 
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Note also 
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that ([5| coincides with the power law tail conjectured by 
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Fyodorov and Sommers [16], Pn{t) 
Coulomb gas.- To derive our main results ( |2|3|4|5[ ), we 



use the Coulomb gas method, originally introduced by 
Dyson [57] . Recently, this method has been suitably 
adopted and successfully used in a number of differ- 
ent contexts : e.g. the distribution of the conduc- 
tance of chaotic cavities [2BH3D], or the quantum en- 
tanglement in a random bipartite state |31H33| . Our 
starting point is to rewrite the joint distribution ([!]) 
of the rescaled rates X{ — "fi/N as a Gibbs measure, 
P(li, • • • ,7jv) oc exp-i/37V 2 £ [p], with the "energy" 8{p] 
expressed as a functional of the density of the rescaled 
rates p(x) — i 2i=i $( x ~ x i)- The energy reads 



£ [p] = / dx (x — hxx) p(x) 



dxdx' p(x) p(x') In \x — x | (7) 



The rescaled time-delay is s — Ar-yv = J2i l^ 1 (i- e - the 
DoS of the cavity in appropriate units s = v(E)A). In 
the limit N —> oo, the density p{x) may be treated as con- 
tinuous and the distribution Pn{t = s/N) can be derived 
via a saddle point method. The optimal (saddle point) 



distribution minimizes Q with two constraints: normal- 
isation J dx p(x) = 1 and J ^ p{x) — s. This requires to 
minimize the "free energy" J-[p] = £{p] + po{ f dx p(x) — 
l) + jiti ( J ^ p(x) — s) , where po an d p\ are two Lagrange 
multipliers that enforce the two constraints (we neglect 
the subdominant contribution of entropy [34). Setting 
the functional derivative x 5 ! F . = gives 



Po + x o ~ ha x + 
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dx p(x) In | x — x a I = . (8) 



where we assume that the optimal density has support 
over the interval xq € [a, b]. Deriving once more with 
respect to xq gives 
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where ~f. represents the principal part. This equation ex- 
presses the force balance at equilibrium, for any charge 
at xq € [a,b], between the confining force — V^(x) com- 
ing from the effective potential V c s{x) — x — lnx + 
and the Coulomb repulsion force from other charges. We 
denote by p*(x;s) the solution of ([9]). The time-delay 
distribution then takes the scaling form 
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where the large deviation function is $-(s) = 
£[p*(x;s)] — £ [p* (x; 1)] (note that when the two con- 
straints are fulfilled, J-[p*\ — £ [p*])- The term £[p*(x; 1)] 
emerges from the normalisation of ([T]), obtained by solv- 
ing the same equation in the absence of the second con- 
straint, i.e. for pi = 0, which we will show to coincide 
with s = 1. Using ([8]), we may rewrite the energy of the 
optimal distribution as 
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x — In x + Xq — In xo 



dxp*(x;s) (11) 
In \x — xq\ 



Optimal distribution. - The integral equation (|9| may be 
solved thanks to a theorem due to Tricomi [35] . We find 
the optimal distribution 
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where the three parameters a, b and c = p\j\fab can 
be found by solving the three algebraic equations ob- 
tained by imposing the vanishing of the density at the two 
boundaries and the condition J ~p*(x;s) — s. These 
equations are conveniently written in terms of the vari- 
ables v — \fab~ and u = \Jajb. A few steps of algebra 
shows that u solves 



a(u) = (1— u) : 



16u 3 + 2u 2 + I6u - 1) 



16m 2 (3m 2 - 2m + 3) 



(13) 
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Then v, /xi and c are given by v — 2u 3 " 1 _^ 2 "t 3 , Ml — 



-4n 2 ( " - 6 "+ 1)(3 " 4 " 2M+3) and c = ^ 

(1 — U J ) 4 D 
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Mosi probable values.- We first analyse the distribution 
P/v(t) in the vicinity of its maximum. £[p* (x; s)] is min- 
imised, i.e. Pn( t ) is maximized, by removing the con- 
straint J Q b — p(x) — s, i.e. by setting /xj = 0. For con- 
venience we introduce the roots x± = 3 ± 2y2 of the 
polynomial u 2 — 6u + 1. For /xi = 0, Eq. ( 13 1 has solution 
u = \J x-/x + = X- with v = 1 and s = 1 and conse- 



quently a 



= 0.171... and 6 = 



5.828... In this 



case we recover the Marcenko-Pastur (MP) law 
1 
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Expansion of Eq. ( 13 ) around the MP point leads to s — 
1 ~ — ^=(u — X-), hence v ~ 1 + f%(u — £-) — 1 — 
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i(s — 1). The corresponding 



energy (111 may be conveniently obtained by choosing 
xq = 1 : we see that the first term is quadratic j(s — l) 2 ; 
we check numerically that the remaining integral term is 
constant, equal to £[p ir (x\ 1)] = 3 — 2 In 2, up to higher 
order corrections [numerical fit gives a correction \{s — 
l) 3 ]. Therefore we conclude that $-(s) — \{s — l) 2 , 

s~l 

i.e. Eq. ([3| (the parabolic behaviour is compared to the 
numerical calculation of the integral ( 11 ) in Fig. [31. 
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FIG. 2: (color online) The optimal density of eigenvalues for 
different values of s ; when s increases, the density eventually 
freezes to the MP law (dashed line). 



Large deviations for s - 
gives u = 1 — \/2s + s 







0.- Expansion of ( 13 ) for s 
0(s 3 / 2 ), hence v = \ + O(s ). 
The support of the distribution is given by a = - [l — 
^/Ts + 0(s)] and b = 1 [l + + 0(s)] (the Lagrange 
multiplier is /xi = +C'(s~ 1 )). The optimal distribution 
ressembles the semi-circle law centered around l/s : 
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This was expected : when s — > 0, the eigenvalues {xt} of 
the Wishart matrix are constrained to be very large and 



they do not feel the spectrum boundary at x = 0. Hence, 
their distribution coincides with the Wigner semi-circle 
law for the usual Gaussian ensembles of random matrices. 
The energy may be conveniently calculated by choosing 
xq = 1/s ; this makes the first term of (111 vanish. The 



leading order of the integral term is straightforwardly 
calculated from (151 : we deduce $_(s) ~ r + I Ins — 
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|(1 — In 2), thus proving pi). The factor exp— ^ is m 
perfect agreement with the exact results for N = 1 & 2 
mentioned earlier. 

Large deviations for s ^ 1 - Freezing transition.- As 
s increases, it eventually reaches a finite value corre- 
sponding to the maximum of the function a(u) (inset 
of Fig. |j), at u c = § [1 + 2(2 1 / 3 - 2 2 / 3 )] = 0.115... giv- 
ing s c = a(u c ) = ^ttJ^l^) = 1-1738.... Then 
a = — c, which leads to a somewhat unusual form 
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For s > s c , (13 1 has no longer physical (real) solutions. 
In this case, the saddle point turns out to have a differ- 
ent solution where a single isolated charge, say at Xi, 
splits off the main body of the density and carries a 
macroscopic weight (see Fig. [TJ. A similar scenario oc- 
curs in the study of quantum entanglement in random 
bipartite state [3"TM33j . We decompose the density as 
p(x) = ±S(x-Xi)+p(x) where p(x) = ^ Yh>\ s (x-x l ) 
is still treated as a continuous density. The energy 



£[ P ]=£[p] 



X\ — lnxi 
N 
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N 



dx p(x) ln(x — xi) (17) 



must be minimized under the two constraints J dx p(x) = 
1 — 4? and f dx = s — -J^- . This leads to the two 

JV J x Nxi 

equilibrium conditions 
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V Xq € [a, b] and X\ < a. We show that a consistent pic- 
ture is the freezing of the density p(x) while the isolated 
charge goes to zero x% —> 0. When N —> 00, the r.h.s. of 
( 19 1 reaches a constant value as x\ — > ; so does the l.h.s. 
iff pi ~ —x\ — > _ . Hence the solution of (18) is the MP 
law: p*(x; s) = p*(x; 1)+0(N~ 1 ). The rescaled time de- 
lay splits into the contribution of the isolated charge and 
of p as s = -jyt- + lj i.e. x\ — l/[N(s — 1)]. In fact this 
analysis holds for any s > 1 (and not only s ^ s c ) : the 
energy (17) of this new phase coincides with the energy 



of the MP solution, up to 1/N corrections. Therefore for 
1 < s ^ s c we have found another phase with a lower 
energy, which shows that the branch obtained previously 



(with compact solution ( 12 ) over [a, b] for s < s c as well 
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as ( p~6[ ) for s = s c ) is actually metastable (Fig. [31). In 
the (thermodynamic) limit N — > oo, the energy of the 
gas vanishes for all s > 1, while for s < 1, it behaves as 
| (1 — s) 2 as mentioned earlier (Fig. [3]). This then results 
in a second order phase transition at s = 1. We call this 
a freezing transition, because for s > 1, energy freezes to 
the value in the thermodynamic limit and also the bulk 
density freezes to the MP distribution. 

One can analyse more precisely this new frozen phase 
by computing the 1/N corrections to the energy. For 



large enough s, Eq. (17) is dominated by the logarithmic 

filn[iV( S - 



term 



jvlnxi, i.e. £[p*{x;s)\ ~ (•••) + £ln 

We get the power law tail Pn{t) ~ (s— 1) ~' , where 
is some exponent of order N° introduced in order to 
account for N~ 2 corrections to £[p]. This exponent may 
be determined as follows : when tw > 1/N, most of the 
proper times are described by the frozen density (the MP 
law), i.e. n G [x-/N, x + /N] for i > 1 with J2i>i 1/ T i = 
1, while one proper time becomes much larger and carries 
a "macroscopic" contribution, t\ = s — 1 = iWw — 1. In 
the scattering problem, this is interpreted as the large 
contribution of a narrow resonance. Writing Pjv(t) = 
J dji ■ ■ -d<y N 5(Nt - I/71 - l)P(7x,--- ,7jv) and using 
leads to = 2, hence Eq. 

A more precise analysis of Eqs. ( 18|19 1 leads 
to introduce the large deviation function 3>+(s) = 
N (£ [p*(x; s)]— £ [p*(x; 1)]) — In N giving the scaling form 



Pn(t) 



N 2 ex p 



■^ + (Nr) forr>^ (20) 



One obtains that (s) = ln(s -l)-l-21n2 (c.f. 

inset of Fig. |3| . The local minimum at s = 2 is related 
to ( |4| while the logarithmic behaviour to the power law 
tail ([5]). For finite N, the energy functions characterising 
the two phases cross for s — sn such that $_(sat) = 
^[$+(sat) + IniV]. Using the limiting behaviours for 
s — > 1, we obtain the finite N correction to the position 
of the phase transition : sjy — 1 + (4/iV) 1 / 3 . 

Conclusion.- In summary, by using a Coulomb gas ap- 
proach, we have analysed the large deviation functions 
controlling the Wigner time-delay distribution in the 
limit of large number of conducting channels. We have 
shown that the distribution exhibits a rich structure. In 
particular, its power law tail is related to a freezing tran- 
sition in the Coulomb gas, corresponding to large contri- 
butions to tw of resonant states in the original scattering 
problem. We have also performed a Monte-Carlo simu- 
lation of the Coulomb gas up to 1600 charges and found 
good agreement with our analytical results (details will 
be published elsewhere). 

Several questions remain open : (i) a more precise 
treatment of 1 /N corrections would be desirable, (it) The 
starting point of our calculation, Eq. 0, describes the 
usual random matrix ensembles ; the distribution of t\v 
was also obtained in [37] for a chiral-GUE ensemble when 
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FIG. 3: (color online). Large deviation function $-(s) (i.e. 
rescaled energy of the gas). The freezing transition takes 
place at Soo = 1. The metastable branch terminates at 
s c = 1.1738... Inset : Large deviation function $+(s) [i.e. 
1/N correction to the rescaled energy]. 



N = 1. Extension of our analysis to such cases would be 
certainly interesting, in particular with the growing in- 
terest in the study of new symmetry classes of disordered 
systems. 
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